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This paper clarifies the microscopic nature of the staggered scalar order, which is specific to even number 
of / electrons per site. In such systems, crystalline electric field (CEF) can make a singlet ground state. As 
exchange interaction with conduction electrons increases, the CEF singlet at each site gives way to Kondo 
singlets. The collective Kondo singlets are identified with itinerant states that form energy bands. Near the 
boundary of itinerant and localized states, a new type of electronic order appears with staggered Kondo 
and CEF singlets. We present a phenomenological three-state model that qualitatively reproduces the char- 
acteristic phase diagram, which have been obtained numerically with use of the continuous-time quantum 
Monte Carlo combined with the dynamical mean-field theory. The scalar order observed in PrFe4Pi2 is 
ascribed to this staggered order accompanying charge density wave (CDW) of conduction electrons. Ac- 
curate photoemission and tunneling spectroscopy should be able to probe sharp peaks below and above 
the Fermi level in the ordered phase. 

KEYWORDS: Kondo effect, crystalline electric field singlet, charge density wave, dynamical mean-field 
theory, quantum Monte Carlo, PrFe4Pi2 



1. Introduction 

In interacting electron systems with two /-electrons per site 
as in some Pr and U compounds, the distinction between itin- 
erant and localized characters of electrons is not as explicit 
as in Ce or Yb compounds with /' or /'^ configuration, re- 
spectively. The difference comes from the Kramers degener- 
acy; the ground state of each non-Kramers ion such as Pr^^ or 
U^^ can be a singlet even without Kondo effect provided the 
CEF favors the non-degenerate ground state. On the contrary, 
the Kramers degeneracy in /' or /'^ can only be removed 
by exchange (Kondo) interaction with conduction electrons. 
The dichotomy between itinerant and localized characters of 
/-electrons leads to intriguing phenomena in heavy fermion 
systems. In particular, mysterious ordering phenomena have 
been found in systems such as PrFe4Pi2 and URu2Si2.''^' In 
the case of PrFe4Pi2, the order keeps the point-group symme- 
try around each Pr, and is called the scalar order. 

A typical situation in /^ systems is the case where the 
ground state is the CEF singlet. If some excited CEF levels are 
involved by interaction with conduction electrons, the ground 
state becomes non-trivial. Let us consider the case where CEF 
states form a quasi-quartet composed of a singlet ground state 
and a first-excited triplet. We represent the singlet-triplet lev- 
els at each site / in terms of two pseudo spins'"'^' Syi with 
y - 1,2. Correspondingly, we introduce two bands for con- 
duction electrons that interact with pseudo spins with the 
same orbital index y by the exchange interaction / > 0. If 
/ electrons interact strongly with conduction electrons, the 
ground state can be a collective Kondo singlet involving CEF 
triplet at each site. In this case, / electrons acquire itinerancy 
because of the Kondo effect. Competition between these two 
singlets, in other words between itinerant and localized char- 
acters of /-electrons, gives rise to rich physics with exotic 
ordered phases. 

The same kind of competition has in fact appeared already 
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in the two-impurity Kondo problem**' where two / electrons 
are spatially separated. In our case the two / electrons are 
on the same atomic site. There are also many theoretical dis- 
cussions in /- impurity systems with the CEF singlet.''"'-'' In 
considering possible ordering in non-Kramers lattice, we can 
draw analogy from the two-impurity Kondo problem and the 
/^ impurity problem. 

In this paper, we provide intuitive picture how the com- 
petition between CEF and Kondo singlets leads to the new 
electronic order that has been proposed recently.'"*' This 
novel electronic order has originally been studied by elabo- 
rate numerical work that combines the continuous-time quan- 
tum Monte Carlo method (CT-QMC)''' and the dynamical 
mean-field theory (DMFT),'**' which is extended for two sub- 
lattices corresponding to staggered order '^•^"' In addition to 
reviewing some representative physical quantities that char- 
acterize the staggered order, we introduce in this paper a toy 
model with three states per site, and clarify the nature of the 
phase diagram, including the tricritical point. We also sug- 
gest a new experiment that can test whether our microscopic 
picture applies to the actual electronic order in PrFe4Pi2. Ac- 
cording to our scenario, each conduction band shows a pair of 
sharp features that develop in the staggered singlet order. 

2. Analogy to two-impurity Kondo problem 

In order to describe the characteristics of our model, we 
first draw analogy to the two-impurity Kondo problem. There 
are two extreme cases concerning the magnitude of the in- 
tersite exchange / and the Kondo exchange J both of which 
are taken to be positive (antiferromagnetic): (i) In the case of 
I » J, the impurity spin pair form the singlet, and J acts as 
weak perturbation; (ii) In the opposite case / <sc 7, each im- 
purity spin forms the Kondo singlet by trapping a conduction 
electron. Then / acts as weak perturbation. In the compet- 
ing case I ~ J, the energy difference between the two kinds 
of singlets becomes smaller than either of I or J. The differ- 
ence becomes simply |7 - /| if the width of conduction band 
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is negligible. It is known from detailed analysis^ that the two 
different singlets have the anti-crossing around J ~ I in gen- 
eral, but they become degenerate if the particle-hole symme- 
try is present in the conduction band. In the former case, as the 
Kondo exchange increases, the localized character of / elec- 
trons, represented by the pair singlet, changes continuously to 
itinerant character, represented by two Kondo singlets. 

We regard the /- configuration in Pr^^ or U'*^ ion as the 
short-distance limit of two Kondo centers. Instead of two real 
spins that tend to form the triplet by the strong intra-atomic 
exchange, we consider two pseudo spins that have antiferro- 
magnetic exchange interaction / (> 0). Then the triplet level 
is located higher by / = A relative to the singlet. The relevant 
Pr and U compounds are modeled as the periodic lattice of 
these pseudo spins. 

3. Generalized Kondo lattice with two conduction bands 

Although the actual Pr and U compounds have compli- 
cated band structures and interactions, we consider the sim- 
plest model that contains the essential feature of the compe- 
tition between itinerant and localized characters. We assume 
the presence of two conduction bands indexed by y (= 1,2) 
each of which interacts with a localized pseudo spin Sy, at site 
/. The model is given by 

kycr iy 

+ A_^Si,-S2/, (1) 

where s^yt denotes the spin of conduction electrons at site /. 
The CEF splitting A is simulated by the coupling between 
pseudo spins Si, and 52/ as shown in the third term.^'^' This 
model at half filling has been investigated in one dimension 
by the density-matrix renormalization group. ^'^ In this paper 
we concentrate on the case where each site has one conduc- 
tion electron on the average. This means that each conduction 
band is quarter-filled if the bands y - 1 , 2 are equivalent. 

We use the dynamical mean-field theory (DMFT) for 
setting up an effective impurity model, and employ the 
continuous-time quantum Monte Carlo (CT-QMC) to solve 
the impurity model. In numerical calculation, the density of 
states of the conduction band is taken to be 



Pc{s) = 




where the origin of energy is at the center of the band. The 
density of states corresponds to the tight-binding model in 
the hypercubic lattice in infinite dimensions. Unless other- 
wise stated, the parameters to be used in the following nu- 
merical results are taken to be: D - 1,7 = 0.8, A - 
0.2. The Kondo temperature in the single impurity Kondo 
model becomes Tk ~ 0.1 according to the formula Tk = 
Dexp{-l/[yp,Ou))]) 

Since the method of calculation has been explained in detail 
in our previous papers, we do not repeat the explanation 
here. Instead, we describe qualitatively the ground state in the 
strong-coupling limit. Figure 1 illustrates the spin and charge 
densities in the ordered phase. As in the two-impurity Kondo 
problem, the competition between J and A in eq.(l) gives rise 
to a small energy scale |/ - A| in the strong-coupling limit. 



Kondo singlets 
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1# # ^ 

CEF singlet 

Fig. 1. (Color online) Illustration of the ordering pattern in the strong cou- 
pling limit. The thick arrows represent localized spins, while objects on the 
two lines represent conduction electrons with spin. 

If the attractive intersite interaction between these two kinds 
of singlets exceeds the difference \ J - A|, the staggered order 
becomes energetically favorable. By precise numerical work 
that fully includes the kinetic energy of conduction electrons, 
we have confirmed that this staggered order is indeed realized 
with certain set of parameters. In the following we describe 
the salient features of this new kind of order. 

4. Phase diagram 

In our calculation the whole lattice is divided into sublat- 
tices A and B. The most fundamental quantity is the site- 
diagonal component G^^(ie„) of Matsubara Green function 
with additional sublattice index /I = A, B. We first derive the 
Matsubara Green function using CT-QMC plus DMFT, and 
then other quantities such as the density of states by analytic 
continuation to real frequencies. If there is no electronic or- 
der, local Green functions do not depend on sublattices. The 
phase transition is the point where the Green function for the 
A sublattice begins to be different from that for B. Another 
important quantity is the two-particle Green function of con- 
duction electrons, and that of / electrons representing local- 
ized spins. From the two-particle Green function with imagi- 
nary frequencies, we obtain the dynamical susceptibilities by 
analytic continuation to real frequencies. On the other hand, 
spin correlation functions are obtained by direct evaluation of 
the statistical average. 

A second-order transition is signaled by divergence of a 
susceptibility that corresponds to the response of the order pa- 
rameter In our case, an element of the staggered charge sus- 
ceptibility which has symmetric orbital combination di- 
verges at the highest T . We represent this element as x^dQ)- 
Figure 2 shows the phase diagram obtained from the diver- 
gence Q^x^ciO)- There appears a region where the first-order 
transition occurs. This is detected from the hysteresis in phys- 
ical quantities as we vary T or A. For example, as A is de- 
creased from A - 0.5 at T = 0.015, 1/a'+c(6) tends to zero 
around A = 0.41. On the other hand, starting from the or- 
dered phase, 1/a-+c(6) is finite at T = 0.015 and A = 0.41, 
and decreases to zero around A = 0.44. Hence, the interven- 
ing region 0.41 < A < 0.44 corresponds to metastable state at 
r = 0.015. The metastable region shrinks with increasing T , 
and terminates at the tricritical point at T ~ 0.2. Because of 
the numerical difficulty, the phase boundary around A ~ 0.46 
has been obtained only crudely. 

We find that the homogeneous susceptibility for the quan- 
tity 

X X "^"^rWcyfr, (3) 
<T=±\ y 
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Fig. 2. Phase diagram with one conduction electron per site. Here and in all 
numerical results, the exchange interaction is set to 7 = 0.8. The tricritical 
point separates the second- and first-order transitions. Estimated upper and 
lower metastability limits in the first-order transition are indicated. See text 
for the phase labeled "Magnetic". 



with T\ - \ and to - -\ diverges at T < 0.01 inside the 
scalar-ordered phase. The order parameter breaks the time re- 
versal even without finite magnetization. Hence this order pa- 
rameter is interpreted as a kind of magnetic octupole. Detailed 
study of this order will be reported elsewhere. 

Let us introduce a phenomenological model that is con- 
venient to interpret the phase diagram. Taking the strong- 
coupling limit, we consider the following three states for each 
site /: (i) the Kondo state with two conduction electrons with 
energy £2; (ii) the CEF singlet state with no conduction elec- 
trons with energy £0; (iii) the CEF singlet state with one con- 
duction electrons with energy £1, which is taken to be 0. The 
states (i) and (ii) are nondegenerate, while the state (iii) is 
four-fold degenerate since the conduction electron has spin 
and orbital degrees of freedom. 

The variable cr, is assigned for each state with a value: (i) 
(T; - 1 ; (ii) o", = - 1 ; (iii) cr, = 0. Furthermore, the constraint 



N 



0, 



(4) 



/■=! 



is imposed to describe the situation of one electron per site 
on the average. This constraint is accomplished by adding the 
Lagrange multiplier term hYjiCi to the Hamiltonian. In the 
mean-field theory, we obtain the partition function and the 
average {cr}^ for each sublattice A -A, B. Assuming (cr)^ s 
cta > and {o-}^ s ctb < 0, we obtain for A =A, 



Za = exp [-/3(E2 +h + /Tctb)] + 4 
+ exp [-/3(£o -h- Ka-^)] , 



(5) 



where (> 0) is the interaction between the nearest neighbor 
sites, and is given in terms of the hopping f as A' ~ 0{zt'/6E) 
where z is the number of nearest neighbors. The excitation en- 
ergy for the state (iii) is written as 6E, which is of the order of 
J or A. Note that K describes the effective attraction between 
the states (i) and (ii). We obtain similar result for Zb. 

We parameterize as £2 - E + d and £0 - E-d where 2d is 
estimated as A - / in the strong-coupling limit. The constraint 



of eq.(4) requires h - -d, and eq.(5) is simplified to 
Za = 2 exp (-yS£) co?.\\(j3Kcr) + 4, 



(6) 



where we have used the constraint cta - -ctb = cr. The parti- 
tion functions of both sublattices become equal to each other: 
Zj^ - = Z. Let us consider some limiting cases, which 
clarify the meaning of £. In the limit £ -00, our model is 
reduced to the antiferromagnetic Ising model. This is under- 
stood as exclusion of the state (iii) with «c = 1 ■ In the oppo- 
site limit of E IK ^ 1, all sites become (iii), and no staggered 
order occurs. Hence, the value of EjK controls the relative 
weight of the state (iii), and subsequently the phase transition. 
The free energy F per site is given by 

F^-TlnZ.\K.\ (7) 

The equation of state is obtained from the stationary condition 

dFjda- = as 



(T = - exp(-;8£) sinhO0/:cr) = g{(r). 



(8) 



Nontrivial solution cr is possible only if E < K. 

In the second-order transition, the nontrivial solution be- 
comes infinitesimal at T^, which tends to Tc - K in the limit 
£ — > -00. As £ becomes larger, the state (iii) acquires a finite 
weight. The transition temperature is then given by 

Tc = , (9) 

1 +2exp(yS,£)' 

which is a transcendental equation for Tc - 1 IPc- The resul- 
tant Tc is smaller than K. 

The second-order transition beginning from £ from -00 en- 
counters a tricritical point at £ = £3 with transition tempera- 
ture T - Tj,. For larger £, the transition becomes of first order 
In order to derive the parameters corresponding to this tricrit- 
ical point, we expand the RHS of eq.(8) up to the third order 
in cr. The coefficient of the third-order term becomes zero at 
the tricritical point. We derive the parameters as 



£3 = 0, 7^3 = K/3. 



(10) 



The first-order transition eventually disappears for larger £ 
at £ = Ec- We obtain E^ = A' by considering the low- 
temperature limit of eq.(8). For £ > E^ the model has no 
phase transition. 

Figure 3 shows the phase diagram of the three-state model 
obtained numerically. The temperature T„ giving the upper 
limit of the first-order transition is determined by the condi- 
tion 



g'i(r) = 1, 



(11) 



for cr 9i 0. The lower limit Ti is given by the same equation but 
with cr = 0. As shown in Fig. 3, the lower metastability limit 
appears mostly as a reentrant transition. The thermodynamic 
transition line corresponds to the condition 



F(cr) = £(0), 



(12) 



where cr (^t 0) is the solution of eq.(8). Since we have cr = 1 
with r — > 0, eq.(12) is satisfied only if 



-{K -E) + K/2 = 0, 
which gives E/K - 0.5 as shown in Fig. 3. 



(13) 
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Fig. 3. (Color online) Phase diagram of the three-state model, together with 
illustration of the free energy as a function of the order parameter cr. As 
E increases, the second-order transition changes to first-order one beyond 
the tricritical point at (£, T) = (0, A'/3). The limits of metastability appear 
on both sides (7",, and Ti) of the thermodynamic transition line. 
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Fig. 4. (Color onhne) Average occupation number of conduction electrons 
at each sublattice A or B. 

These results are useful to understand the numerical results 
in Fig.2, which are not in the strong-coupling limit. We in- 
terpret Tk as playing the role of J. Since increase of A > 
causes decrease of Tk, the transition temperature increases up 
to the point where Tk becomes comparable to A. In the phe- 
nomenological three-state model, increase of 2d - A - Tk, 
with E() - E - d kept fixed, means increase of E. Namely, the 
increase of A in Fig.2 corresponds to increase of E. Thus the 
change of second-order transition to first-order one, and even- 
tual disappearance with increasing A can be interpreted by the 
three-state model. It is remarkable that the salient feature of 
Fig.2 is reproduced quahtatively by this simple model. 

5. Static quantities in ordered phase 

5.7 Staggered charge density 

Figure 4 shows the number of conduction electrons per site. 
This quantity has been obtained for each sublattice /I = A or 








T 



Fig. 5. (Color onhne) Local congelation functions (a) between conduction 
and localized spins representing the Kondo effect, and (b) between local- 
ized spins representing the stability of the CEF singlet relative to the triplet. 
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We have fixed the average number - I distributed over the 
two equivalent conduction bands, as seen in Fig.4 above the 
transition temperature Tc ~ 0.032. As temperature T becomes 
lower than Tc, the disproportion occurs between the sublat- 
tices A and B. Their sum «cA + "cB remains constant (= 2) as 
in r > Tf. by construction. Hence the electronic order we have 
obtained corresponds to a CDW state of conduction electrons. 
The CDW is related to the difference of local spin correlations 
as shown next. 

5.2 Local spin-correlation functions 

Let us examine equal-time spin correlations that clarify 
properties of each sublattice with the finite CEF splitting. We 
first consider (Si • Sd), which is equal to (S2 ■ Sci) under the 
present condition. Figure 5(a) shows the temperature depen- 
dence of {Si -Sci), which is enhanced on A-sublattice and sup- 
pressed on B-sublattice. Bearing «cA > "cB in mind, we con- 
clude that the localized spin on A-sublattice forms the Kondo 
singlet. On the other hand, the correlation (Si ■ S2) between 
localized spins is enhanced at B site as shown in Fig. 5(b). 
Hence, the B-sublattice corresponds to the CEF singlet. 

The magnitude of the effective CEF splitting for the CEF- 
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singlet site can be estimated as follows: Taking the effective 
Hamiltonian for the localized states as - ASi -82, we ob- 
tain the susceptibility as Xu ~ ~ 1/(2 A) for the ground state. 
Here A is the effective CEF splitting, which can be estimated 
from the local susceptibility in B-sublattice. The numerical 
result is A = 0.098, which is very close to the original CEF 
splitting A = 0.1. Hence, spatially extended Kondo singlets 
do not significantly affect the magnitude of the CEF splitting. 
Besides, in Figure 5(b), (Si ■ S2) in the low-temperature limit 
is not far from -0.75 expected for the isolated singlet. Hence 
the CEF singlet interacts only weakly with conduction elec- 
trons. Note, however, that the correspondence between A and 
A does not hold for A < 0.05. This is because the CEF-singlet 
site tends to be magnetically polarized near A = 0. 

In this way, this electronic order turns out to be a staggered 
order with the Kondo and CEF singlets, as illustrated in Fig. 1 . 
Except for the strong coupling limit, the number of conduc- 
tion electrons at the CEF singlet site is not zero because the 
Kondo singlets are spatially extended. This fact becomes im- 
portant in understanding the single-particle density of states 
to be presented in the next section. 

6. Dynamical quantities 

6.1 Dynmnical susceptibility 

In our theoretical framework, it is in principle possible 
but numerically tedious to derive the momentum-dependent 
dynamical susceptibility. Hence, prior to full derivation of 
the dynamical susceptibility in the future, this paper presents 
some results only for the local dynamical susceptibility that 
corresponds to the momentum average. In skutterudites such 
as PrFe4Pi2, the triplet wave functions IF,) are mixture of the 
cubic r4 and Vs states because of the lower local symmetry 
Th. Namely the wave function is given by 



60 



IF,) = V^|r4> + Vr^iFs). 



(15) 



Then matrix elements of the magnetic dipole between singlet 
and triplet levels are fixed for given w. Note that the limit 
w = makes the dipole moment vanish. 

Figure 6 shows the results related to the dynamic mag- 
netic susceptibility ;t'y(w) above and below the transition tem- 
perature. Here we have used the parameter w - 0.64 for 
the numerical calculation. At T = 0.04, which is above the 
transition temperature Tc ~ 0.032, the Lorentzian spectrum 
emerges characterized by the quasi-elastic components. As T 
decreases below T^-, the first inelastic feature appears with the 
peak at coi ~ 0.14. With further decrease of T, the second 
inelastic peak also develops at llio ~ 0.18. By looking at sub- 
lattice components, we have identified the first inelastic peak 
at a)[ comes from the CEF singlet sites, while the second peak 
at ai2 from the Kondo singlet sites. 

6.2 Single -particle spectrum 

The momentum average of single-particle dynamics can be 
seen in the density of states. Figure 7 shows the renormalized 
density of states given by 



,1=A,B ycr 



-{s + i(5). 



(16) 







X 



T = 0.040(para) 
T = 0.025 
T = 0.010 




Fig. 6. (Color online) Local dynamical magnetic susceptibility at different 
temperatures. Below the transition temperature ~ 0.032, A and B sub- 
lattices have different spectra, as shown by thin lines. The Kondo-singlet 
site has the broader and higher-energy peak, while the CEF singlet site has 
the shaper and lower-energy peak as most clearly seen at T = 0.010. 
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for J - 0.8 and A = 0.2. Here the local Green function 
Gy^ie+i6) of conduction electrons is derived by analytic con- 



Fig. 7. (Color online) Density of states for conduction electrons at different 
temperatures. The Gaussian density of states of non-interacting case ( J = 
0) is also shown. The sharp peak below the Fermi level at 7" = 0.015 < Tc 
comes from the Kondo-singlet site, while the other peak above the Fermi 
level from the CEF singlet site. 



tinuation of Matsubara frequencies. Under the present condi- 
tion, the Green functions does not depend on the labels y and 

O". 

In the disordered state {T - 0.04), the density of states 
is finite at the Fermi level, indicating the metallic behavior. 
The large dip structure around cj ~ 0.3 is interpreted as a 
kind of hybridization between conduction electrons and the 
localized states by the Kondo effect. Although there is no real 
hybridization because the pseudo spins do not have charge 
degrees of freedom, strong renormalization by the Kondo ef- 
fect gives rise to an electronic state that allows this interpre- 
tation.-"*' At r = 0.015 which is a temperature inside the or- 
dered phase, the double peaked structure arises as in Kondo 
insulators. In the present case, however, each peak is asso- 
ciated with different sublattices. Namely, the sharp peak be- 
low the gap comes from the sublattice for the Kondo singlet, 
while the peak above the gap is due to the CEF singlet site. 
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Fig. 8. (Color online) Renormalized energy bands schematically shown in 
one dimension. A pair of original conduction bands hybridize with Kondo 
states, and end up with six energy bands in the reduced Brillouin zone in 
the ordered phase. We have used e = - cos k + const, for the spectrum of 
bare conduction bands. The flat parts in the renormalized bands appear as 
peaks in p(e). With one conduction hole per site, the Kondo-singlet site has 
a peak above the Fermi level, while the CEF singlet site has another peak 
below the Fermi level. This con'esponds to the reversal of energy from the 
result shown in Fig. 7. Furthermore, the two conduction bands here are not 
equivalent, as indicated by different colors. 

Hence, this double peaked structure clearly shows the differ- 
ence of the occupation number between the Kondo and CEF 
singlet sites. We emphasize, however, that the Kondo effect is 
responsible for both peaks, as in the case of Kondo insulators. 

The origin of the insulating behavior is explained as fol- 
lows. Each conduction band has one conduction electron per 
doubled unit cell in the ordered state. Provided that the lo- 
calized spin at the Kondo-singlet site participates to the con- 
duction band, each band is filled by two electrons. Then the 
system can be an insulator. In this viewpoint, the staggered 
Kondo-CEF singlet order may be regarded as alternating itin- 
erant and localized sites of / electrons. 

7. Discussion and conclusion 

Let us finally discuss the relevance of the present model 
to real systems such as PrFe4Pi2. We have to consider two 
important effects in PrFe4Pi2: (i) the two conduction bands 
are not equivalent; (ii) there is one conduction hole, rather 
than one conduction electron, per Pr ion. As a result of (ii), 
the CDW gives more holes to the Kondo singlet site than the 
CEF site. Then the Kondo singlet site has the peak above the 
Fermi level for the local density of states of conduction elec- 
trons, and the CEF singlet site has the peak below the Fermi 
level. The ensuing renormalized energy bands are illustrated 
schematically in Fig. 8, where the effect (i) is also taken into 
account. The number of itinerant electrons is 8 (= 6 H- 2) per 
unit supercell in the ordered phase, out of which two electrons 
originate from the Kondo-singlet site with configuration. 
Hence the ground state is a semimetal or semiconductor, de- 
pending on the degree of asymmetry in conduction bands. 

The CDW of conduction electrons should lead to distortion 
of crystals keeping the local point-group symmetry. The dis- 
tortion has indeed been observed by diffraction experiments 
by X rays and neutrons. In inelastic neutron scattering on 



PrFe4Pi2, clear inelastic features appear only below the tran- 
sition temperature.^^- If a small amount of Pr is replaced 
by La, the scalar order disappears, and no inelastic peaks 
are visible even at low enough temperatures.^''* The results 
shown in Fig. 6 reproduce the characteristic features observed 
in PrFe4Pi2, giving further support to relevance of our model. 

The clearest confirmation of the present scenario will be 
provided by accurate spectroscopy such as the photoemission 
and tunneling spectroscopy in PrFe4Pi2. Namely, only below 
the ordering temperature, sharp peaks of electronic density of 
states develop both below and above the Fermi level. Photoe- 
mission probes the sharp feature below the Fermi level. Above 
the Fermi level, on the other hand, the hybridized heavy- 
electron band may be probed by tunneling spectroscopy, 
which has been utilized effectively in URu2Si2.^°' 

In summary, we have described the new staggered singlet 
order for non-Kramers /-electron systems, and discussed its 
salient features in static and dynamic quantities. For the phase 
diagram, we have proposed a simple model and reproduced 
the first- and second-order transitions depending on the rela- 
tive stability of CEF and Kondo singlets. The tricritical point 
measures the strength of the intersite interaction that is re- 
sponsible for the staggered order. The scalar order observed 
in PrFe4Pi2 is ascribed to this singlet order. 

This work has been supported by a Grant-in-Aid for Sci- 
entific Research No. 20340084, and another (No. 20102008, 
Innovative Areas "Heavy Electrons" ) of The MEXT, Japan. 

1 ) For a review, see: Y. Kuramoto, H. Kusunose, and A. Kiss: J. Phys. Soc. 
Jpn. 78 (2009) 072001. 

2) For comparison between PrFe4Pi2 and URu2Si2, see: E. Hassinger et 
al: J. Phys. Soc. Jpn. 77 (2008) Suppl. A, p. 172. 

3) A. Kiss and Y. Kuramoto: J. Phys. Soc. Jpn. 75 (2006) 103704. 

4) O. Sakai et al: J. Phys. Soc. Jpn. 76 (2007) 024710. 

5) J. Kikuchi, M. Takigawa, H. Sugawara, and H. Sato: J. Phys. Soc. Jpn. 
76 (2007) 043705. 

6) R. Shiina: J. Phys. Soc. Jpn. 73 (2004) 2257. 

7) J. Otsuki, H. Kusunose, and Y. Kuramoto: J. Phys. Soc. Jpn. 74 (2005) 
2082. 

8) I. Affleck, A.W.W. Ludwig, and B.A. Jones: Phys. Rev. B52, (1995) 
9528, and references therein. 

9) Y. Kuramoto: in Transport and Thermal Properties of f-Electron sys- 
tems, (G. Oomi, H. Fujii, and T. Fujita eds.. Plenum Publ. Corp., 1993) 
p.237. 

10) Y. Shimizu, O. Sakai, and Y. Kuramoto: Physica B 206-207 (1995) 135. 

1 1) M. Koga and H. Shiba: J. Phys. Soc. Jpn. 65 (1996) 3007. 

12) K. Hattori and K. Miyake: J. Phys. Soc. Jpn. 74 (2005) 2193. 

13) S. Hoshino, J. Otsuki, and Y. Kuramoto: J. Phys. Soc. Jpn. 78 (2009) 
074719. 

14) S. Hoshino, J. Otsuki, and Y. Kuramoto: J. Phys. Soc. Jpn. 79 (2010) 
074720. 

15) A. N. Rubtsov, V. V. Savkin, and A. I. Lichtenstein: Phys. Rev. B 72 
(2005)035122. 

16) R Werner and A. J. Millis: Phys. Rev. B 74 (2006) 155107. 

17) J. Otsuki, H. Kusunose, P. Werner, and Y. Kuramoto: J. Phys. Soc. Jpn. 
76(2007) 114707. 

18) For a review, see: A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- 
berg: Rev. Mod. Phys. 68 (1996) 13. 

19) R. Peters and T. Pruschke: Phys. Rev. B 76 (2007) 245101. 

20) S. Hoshino, J. Otsuki, and Y. Kuramoto: Phys. Rev. B 81 (2010) 1 13108. 

21) S. Watanabe, Y. Kuramoto, T. Nishino, and N. Shibata: J. Phys. Soc. Jpn. 
68(1999) 159. 

22) For a review, see: H. Tsunetsugu, M. Sigrist and K. Ueda: Rev. Mod. 
Phys. 69(1997) 809. 

23) J. Otsuki, H. Kusunose, and Y. Kuramoto: J. Phys. Soc. Jpn. 78 (2009) 



J. Phys. Soc. Jpn. 80 (2011) SAiche-proc 



SAiche-proc-7 



034719. 

24) J. Otsuki, H. Kusunose and Y. Kuramoto: Phys. Rev. Letters 102 (2009) 
017202. 

25) K. Iwasa et ah Physica B312-313 (2002) 834. 

26) K. Iwasa et ah Acta Phys. Pol. B 34 (2003) 1117. 



27) L. Hao et ah. Physica B 359-361 (2005) 871. 

28) J. -G. Park et ah Phys. Rev. B 77 (2008) 085102. 

29) K. Iwasa et ah J. Phys. Soc. Jpn. 77 (2008) 063706. 

30) A.R. Schmidt et ah Nature 465 (2010) 570. 



